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Abstract. A general method is given to derive the current-potential relations in 
anisotropic plasmas. Orbit limit current is assumed. The collector is a conductive 
sphere or an infinite cylinder. Any distribution which is an arbitrary function of the 
velocity vector can be considered as a superposition of many mono-energetic beams 
whose current- potential relations are known. The results for two typical pitch an- 
gle distributions are derived and discussed in detail. The general properties of the 
current potential relations are very similar to that of a Maxwellian plasma except 
for an effective temperature which varies with the angle between the magnetic field 
and the charging surface. The conclusions are meaningful to generalized geometries. 

The Introduction 

The current collection from incoming particle is a fundamental problem in charg- 
ing theory. In this paper the current collection in anisotropic plasmas is studied. 
The analytic expressions for current potential relations which have been used so far 
are derived from Maxwellian distribution. 

Here the basic current potential relation for Maxwellian plasma are repeated 
briefly. These well known results will be compared to that of anisotropic plasma 
frequently in rest part of this paper. The orbit limited current to a spherical con- 
ductor or cylindrical conductor for repelled particles is ( Mott-Smith and Langmuir 
1926; Prokopenko and Laframboise 1977,1980;): 

/ = I D e^ ( 1 ) 

For attractive particle, the charging current is: 

I = Io{l + P) ( 2 ) 

for a spherical conductor and the following: 

/ = / 0 (2/A)(v^+ ( 3 ) 

for a cylinder. 

Charging currents to satellites in anisotropic plasmas have not been discussed 
systematically. The real distributions of the plasmas at synchronous orbit may be 
very different from Maxwellian. For example, the pitch angle distributions in the 
earth magnetic field are very common. In this paper, the current collection in 
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anisotropic plasmas is studied to see how it would be deviate from that in isotropic 
plasma. 

In order to discuss the variety of distributions, all other conditions are assumed 
as simple as possible: 

1) The effect of space charge is neglected. Orbit limited current is assumed. 

2) The distribution function can be an arbitrary function of the velocity vector but 
it does not depend on the spacial location. The currents from two distributions are 
derived in this paper: 
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( 4 ) 


Tj| and T± mean the directions related to the magnetic field, 
b. A source cone or a loss cone in a Maxwellian distribution function. 

3) The probe is a conducting sphere or an infinite cylinder. 

4) VXB force is neglected. The role played by the magnetic field is only to indicate 
the direction of the anisotropy. 

5) Plasma is collisionless. 

The above assumptions are reasonable for satellites at the synchronous orbit. 

A distribution function of plasma and its current to a charged conductor can 
be separated and superposed, the plasma are separated into many monoenergetic 
beams whose charging currents are known. The currents to a sphere conductor 


I p from such a beam are( Mott-Smith and Langmuir 1926 ): 

I p = I 0 { 1 + e<f>/E) = 7T R 2 i{l + e<f>/E) 


( 5 ) 


for attracted particle, the following: 

I p = / 0 (1 - e<f>/E) = 7i" R 2 i{l - e<f>! E) (6) 

for repelled particle with E > e<f> . When E < e<f> , I p = 0 

Here 

g — zL 

P kT 

The ”e”is the magnitude of a electron charge. 
e<f> is the potential energy at the surface, 
k is the Boltzmann’s constant. 

T is the temperature of the plasma 

J 0 is the current to a uncharged probe. I and I 0 are define as total current for the 
sphere and refer to the current per unit length for the cylinder. 

/ 0 = area x . 

m is the mass of a ion. 

E is the kinetic energy of a particle at infinity. 

The i is the current density of the beam. 

R is the radius of the sphere. 
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The equation (4) and (5) simply comes from the energy conservation and the 
angular momentum conservation. 

The basic assumption is that particles which carry maximum angular momen- 
tum and still reach the sphere are the grazing particle. This may not be true for 
attracted particles if the space charge effect is included (Laframboise 1965 Fig 4d). 
Therefore the condition for equation (4) and (5) is that current be orbit limited. 

Similar relations hold for the cylinder: 

I = I Q \J 1 e<j>/E — 2 R + e<£/ E (7) 

for attracted particle and the following: 

I = I 0 \J 1 - e<f>/E = 2 R iy/l-e^E (8) 

for repelled particles with E > e<f> when E < e<f>, I p = 0 

The assumption and derivation of (7) and (8) is the same as equation (5) and (6). 

Caution should be observed: If the beam is not perpendicular to the axis of the 
cylinder, the i in (7) and (8) is not the current density of the beam but is the 
component perpendicular to the axis at infinity. Also the E is not the total kinetic 
energy of the particle but the kinetic energy in the direction perpendicular to the 

axis. Integrating of these beam with weight of distribution function leads to the 

total charging current of incoming plasma. The integral is carried out in the plasma 
frame. The superposition method will be illustrated in more detail in next section. 

The Result of the Current-potential Relations 

It is more convenient to use dimensionless quantities: 

Iq : the current to an uncharged conductor 
I/Io: Dimensionless charging current 
j3 :Nondimensional potential /? = e<f>/kTj_ or e<f>/kT 
k :the Boltzmann’s constant. 

T : the temperature of the plasma 

The current to a attracting sp here from a two temperature plasma. 

The "l” in equation (4) is the current density of the beam at infinity. It is 
equal to the density times the velocity. Now the density should be replaced by the 
number of particles in a infinitesimal volume in velocity space f(v)d 3 (v) . f(v) is 
the distribution function (4) 

The current to the charged sphere is: 

I = f J f f(v)V(l + e<j>/E)d 3 (v)S (9) 
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Here S = n R 2 

d 3 (v) = V 2 d<j> sinO dO dV 

The angle 6 and <)> indicate the directions of particles at infinite. 9 is the polar 
angle of the velocity. 4> is the azimuthal angle of the velocity. All quantities are 
defined in the plasma frame. The distribution function f{v) , the velocity V, and 
the kinetic energy E are the values at infinity (i.e. in the plasma frame). The integral 
is carried out in the plasma frame. 


/ = Nen l iV /„ #/ 0 nnti> 




-imV 2 ( 


2 ( pin 2 9 i c 09*9 \ 


A: Tii 


kT ± 


Define 


B = 


2 A „ ,, 2/ 271\3 


m i 

K(«) 


— Ne R 2 (-)VT{T^ 
m " 


sin 2 0 cos 2 6 

kT]\ + kf[ 
E = \mV 2 


I = B f°° dE f K sin9 dOE {e~ E K{9) + e<f> e~ E K(9) ) 
Jo Jo 


( 10 ) 

( 11 ) 


( 12 ) 


Define I 0 as the current to an uncharged sphere. 


Jo = 2 tt R 2 N\ 


kT\\ 
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Here C = £* 
r li 

D = 2 tg~ l y/C - 1 if T x > Tj[ 

ZY>r„ 

Define the effective Temperature as: 


T, = Tj.( 


yj 1-c 


CD 




(14) 


The current to a charged sphere becomes 

I = J 0 (l + e^/fcT,) (15) 

The dimensionless current are plotted in Fig. 1 as the function of the dimen- 
sionless potential for different values of 2x/T|| which is the indication of the 
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anisotropy. The lower curves correspond to lower Tj| . The curves are simply 
straight lines. The different slopes of the curves is calculated by equation (14). The 
form of current-potential relation is similar to the one of the Maxwellian plasma 
equation (2). The the relation between the current and the potential is linear. 
From equation(5), the current potential relation is linear for a beam to a sphere. 
Therefore the superposition of the beams is also linear. This is clear according to 
equation (13). The change of the charging current is : 

A/ = 7-/ 0 = e^ J f{V)~d 3 v 


The current to a uncharged probe is: 

Jo = / f(V)V d 3 v 


A I/I = e<f> 


i 

f f(V)V d*v 


If we write . 

The change of the current is proportional to the potential. The quantity Tf is 
independent of the potential and has the unit of energy. So we have the definition 
of T f : __ 

f f{V)V 

f f j ^y-J(V)Vd*v 

If the distribution function is Maxwellian. Tf identifies with the temperature. In 
an anisotropic plasma, T f is related to the energy in the direction perpendicular 
to the surface. It is a combination of Tj_ and T\\ . Its values lie between T± 
and Tj| and is more close to T± since Tj_ corresponds to two dimensions 
while Tji only corresponds to one dimension. 


The current to a repelling sphere of two temperature plasma 


Starting from equation (12). I changed the lower limit of the integral and the 
sign before e(f> 


I = B f°° dEE (e~ E K{9) - ecj> <T B KW )sitid d6 
J e<f> 

Here B is defined by equation (10) K (& ) is defines by equation (11) The charging 

current is: 


I = 2n R I 2 
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(17) 


Here G = - j£r , Define /? = e<j>! kT x 
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J 0 is the current to an uncharged sphere defined by equation (12) 

In Fig. 2 the dimensionless current J/Jo is plotted as the function of dimen- 
sionless potential e<j)/kT x Each figure correspond to a different ratio of T x /T\\ . 
The shapes of the curves are very similar to the exponential form of Maxwellian. 
The logarithm scale is used for dimensionless current to show this similarity. 


The current collected by a sphere from a plasma with a loss cone 


For repelled particles, 

J = Ne 7 T f^mdVr x ft* d<f> f^ e " sin6 dO 

9 0 is the angle of loss cone 

J = I om e~$cos9 0 

Here I om is the current to a uncharged sphere in a Maxwellian plasma. 

, / kT j. 

/„„ = / = 4,11^ — 

The Current to an uncharged sphere J 0 is 

Jo = Jom^OS^o ' 

J/J 0 = e -$ (18) 

For attracted Particles 2 

/ = Ne * IS°dV r e^(l + 

r2ir 

/ d(f> sinO dd 

Jo Jo 0 

J = J om (l + ^)cos0 o 
Jo = IomCOS0o 

W* = 1 + if (19) 

The I — <f> relation is exactly the same as Maxwellian plasma. This con- 
clusion can be generalized to any pitch angle distribution if the angle dependence 
is separated from energy dependence, i.e. f(v) = f{9, <f>) x f{E) This result 
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supports the condition that the angle preference of a distribution function does not 
change the I - <f> relation for a spherical conductor if the energy dependence of 
the distribution is the same as the Maxwellian distribution. 


The current to a cylind er in a two temperature plasma. 


Define the coordinates as shown in Fig. 3. The axis of the cylinder is defined 
as the Z direction. X axis is in the plane of B field and axis of the cylinder. The 
magnetic field is parallel to the plane with azimuthal zero. The polar angle of B 
field is 9 o . The azimuthal angle of the velocity is <j>. The polar angle of the velocity 
is 9. || and J_ in The V x is the velocity component in the direction of axis of 

the cylinder. V r is the velocity perpendicular to the axis. 


x e 




j n r°° use for repelled particles, use 0 for attracted particles. 
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Here JO is the zero order first kind of Bessel function E r = 2 m ^r 
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For H > 0 . Here, U = E r /T ± E r is the kinetic energy perpendicular to the 
surface of the cylinder. 


H= 1 + 


sin 2 go(l — C) 


(C - 1) cos 2 9 0 + sin 2 9 0 

T x 


( 22 ) 


C = 


Q 0 is the angle between the B field and the axis of the cylinder 

JS n x e -“l 1+H U0(77u) dU 

7/70 = / 0 °° v /tIe-“( 1+H )/0(Hu) dU 


(23) 
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Here 



in VUT0 use + for attracted particles, use - for repelled particles 

See Fig.4 and Fig. 5 In each figure the curves which change slower with potential 
correspond to higher T\\ . The angle effect is indicated by the H factor in equa- 
tion(22). When the B field is parallel to the the axis of the cylinder, H=0 and the 
IO(Hu) in equation (23) approaches 1. The I — 4> relation is reduced to the 
form of Maxwellian plasma. In this situation the velocity parallel to the axis will 
not be changed by electric field and does not contribute to the charging current at 
all. Therefore Tj| does not appear in equation (23). While 9 0 decreases from 
90° to 0, the motion parallel to B contributes less and less to the charging. Tj| 
becomes less important. 


Current in a Maxwellian plasma within a small source cone 


Iattrac/10 =Ue U Ki(U) 


(24) 


Irepti/Io = U e~ u Ki(U) (25) 

Here U = e <j>/2T sin 2 9o K\ is the Third kind Bessel function of order one. If 9 0 is 
small current-potential relation behaves like a Maxwellian distribution with a lower 
temperature. If plasma come from a direction which is almost perpendicular to the 
surface, The curves behave like a Maxwellian plasma with a higher temperature. 

See Fig.6 and Fig. 7 I/I 0 is plotted as the function of e<f>/kT . A curves 
of Maxwellian distribution is plotted in each figure to be compared with source 
cone. The U factor in equation (25) and (26) shows that 2T sin 2 9o is the effective 
temperature. T e ff ect is less than T when 0 o < 45°. T e ff ec t is greater than T 
when Bo > 45°. The temperature of a Maxwellian plasma can be understood as 
the average of these T e ff ect over all direction. 

I only calculated the I — <f> relation for a flow from a small solid angle, but 
the behavior of a wider source cone can be estimated from the result of the small 
solid angle. If the maximum angle between the axis and a beam within the wide 
solid angle is 60° and the minimum angle is 30° . The I — <j> curve of such 
a source cone should lie between the two curves corresponding to 0° = 30° and 
0 0 = 60° . 


The Discussion 
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Arbitrary Distribution 


There are other distributions such as a double Maxwellian and a monoener- 
getic beam plus a Maxwellian distribution function, for which charging currents are 
very easy to be obtained by using superposition. In general the method in this 
paper applies to any homogeneous distribution function. The only problem left is a 
mathematical one which refers to integrals of equation (5)- (8) in three dimensional 
space of canonical momenta. Sometime integrals have to be carried out numerically. 

The similarity to the Maxwellian plasma 

The general properties of the current-potential curves of pitch angle distribu- 
tions are very similar to that of the Maxwellian plasma The current potential re- 
lations of attractive particle for a sphere is exactly linear as shown in Fig.l . For 
repelled particles the relation are almost exponential for both sphere and cylinder 
as shown in Fig. 2, 5 and 7. 

The importance of the effective temperature 

While the energy of Maxwellian plasma is indicated by the the temperature, the 
energy of a pitch angle distribution in the charging problem refers to an effective 
energy in the direction perpendicular to the charged surface. 

The value of the effective tem perature 

For a pitch angle distribution defined by equation (1), The value of effective 
temperature is an average of Tj. and Tj| . If the surface has no angle preference 
(a conducting sphere) T effect is defined by equation (9). T x contributes more to 
the average since T x indicates the kinetic energy in two dimensions while T| 
is only related to the motion in one dimension. When the surface is parallel to the 
B field, T e ff ee t will deviates from the value of equation (9) and move closer to the 
2j_. When the surface is perpendicular to the B field, T e j fect will approach T\\ 

For the charging of a cylinder from a source cone of a Maxwellian plasma the 
T e ff e ct equals the temperature of the Maxwellian distribution times an angle fac- 
tor. The factor is less than 1 when the source cone makes a small angle with the 
surface. When the source cone is perpendicular to the surface, the T e ff e ct will be 
greater than T. Teffect can exceed 2T. 

Current to a uncharge d surface element 


The current to a uncharged surface has not been studied in section 2. All for- 


82 



mula and Figures are shown in terms of dimensionless currents I / I Q . There is 
not much physics involved in the calculation of I 0 . It is just the random flux 
along the normal direction of the surface. Obviously the current varies with the 
angle between the surface and the preferred direction of the anisotropic distribution 
disregarding the geometry of the whole satellite. The change of the electron current 
and ion current due to I 0 will affect the equilibrium potential. 

The equilibrium potential varies with the orientation of the surface. 


The equilibrium potential varies with the angle between the magnetic field and 
the charging surface. The surface which is parallel to a source cone feel that the 
plasma has less energy; therefore, it would be charged to a lower equilibrium poten- 
tial than the surface perpendicular to the source cone. The equilibrium potential 
depends on the effective temperature. The Tf of a source cone has a upper limit 
which is about 2T. 


The Conclusion 

In case of (4) and source cone, the properties of current-potential relation of a 
conductor in a anisotropic plasma are qualitatively similar to that in Maxwellian 
since the energy distribution is similar. The difference caused by anisotropy is that 
a surface tangential to the preferring direction of anisotropy starts charging with 
a less charging current and reaches a lower equilibrium potential than a surface 
perpendicular to the preferred direction. 
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Fig. 1 The current to a sphere from attracted particles of two temperature 



Fig. 2 The current to a sphere from repelled particles of a two temperature 







